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, ABSTRACT 

■ We use group size haloes, with masses in the range 10 13 < M < 2 x 1O 14 /i _1 M , iden- 

tified with a "friends of friends" (FOF) algorithm in a concordance ACDM GADGET2 
y—{ ■ (dark matter only) simulation to investigate the dependence of halo properties on the 

environment at z = 0. The study is carried out using samples of haloes at different dis- 
tances from their nearest massive cluster halo, considered as such if its mass is larger 
than the upper limit of the above halo mass range (ie., M > 2 x 10 14 /i -1 Mq). We 
find that the fraction of haloes with substructure typically increases in high density re- 
gions. The halo mean axial ratio (c/a) also increases in overdense regions, a fact which 
is true for the whole range of halo mass studied. This can be explained as a reflec- 
[ tion of an earlier halo formation time in high-density regions, which gives haloes more 

""""^ . time to evolve and become more spherical. Moreover, this interpretation is supported 

' by the fact that, at a given halo-cluster distance, haloes with substructure are more 

elongated than their equal mass counterparts with no substructure, reflecting that the 
Q ■ virialization (and thus sphericalization) process is interrupted by merger events. The 

5-^ | velocity dispersion of low mass haloes with strong substructure shows a significant 

c/3 , increase near massive clusters with respect to equal mass haloes with low-levels of 

substructure or with haloes found in low-density environments. The alignment signal 
between the shape and the velocity ellipsoid principal axes decreases going from lower 
•rH , to higher density regions, while such an alignment is stronger for haloes without sub- 

' structure. We also find, in agreement with other studies, a tendency of halo major 

\ axes to be aligned and of minor axes to lie roughly perpendicular with the orientation 

of the filament within which the halo is embedded, an effect which is stronger in the 
proximity of the massive clusters. 

Key words: galaxies: clusters: general - galaxies: haloes - cosmology: dark matter - 
methods: N-body simulations. 
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1 INTRODUCTION 

According to the current Cold Dark Matter (CDM) 
paradigm, haloes emerging from a Gaussian primordial den- 
sity fluctuation field, assemble through gravitational pro- 
cesses to form larger systems which eventually virialize. 
These structures evolve in a hierarchical fashion aggregat- 
ing smaller mass systems, flowing out of voids and along 
filaments, giving rise to deep potential wells, the cluster of 
galaxies. The role that the environment plays in modifying 
the properties of the smaller systems, such as galaxies, is 
being exhaustively studied and it is well known that many 



of the observed galaxy properties correlate strongly with 
environment (eg. Dressier 1980; Goto 2003; see Boselli & 
Gavazzi 2006 for a recent review). The properties of the 
galaxy group-size haloes within which they are embedded 
could also vary as a function of environment, since merg- 
ers and tidal interactions are more probable in high density 
environments. 

The proximity of a galaxy or a group-size halo to a mas- 
sive attractor, like a cluster, and the corresponding strong 
gravitational interactions not only with the cluster itself but 
also with its local surrounding, which is denser near the 
cluster, might affect halo properties such as, among others, 
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shape, size, concentration, orientation, velocity dispersion, 
amount of substructure and internal alignments (for exam- 
ples of such observational evidences, see Schuecker et al. 
2001, Plionis & Basilakos 2002, Plionis 2004, Martinez & 
Muriel 2006). Quantifying such effects in numerical simu- 
lations and understanding their significance could help un- 
derstand the physical processes that act to determine the 
properties of galaxies as a function of their environment. 

However, one should also remember that differences be- 
tween galaxy and group haloes, in high and low density 
regions, could also arise as a natural consequence of cos- 
mological initial conditions, like halo formation time (eg., 
Gottlober, Klypin & Kratsov 2001; Sheth & Tormen 2004) 
or halo spin generation efficiency as a function of local den- 
sity (Lee 2006). 

Various recent studies have applied environment detect- 
ing algorithms in an attempt to characterize the diversity 
of cosmic environments from voids, to walls, filaments and 
clusters and thus facilitate the study of environmental effects 
on galaxy, group and cluster properties. Such algorithms are 
based on a variety of pattern recognition techniques from the 
simplest local overdensity measures to more elaborate tech- 
niques based on second-order local variations of the density 
field (eg. Pimbblet 2005; Colberg et al. 2005; Stoica et al. 
2005; Aragon-Calvo et al. 2006; Hahn et al. 2006) 

Lemson & Kauffman (1999), explored the effect of envi- 
ronment on different halo properties like their mass function, 
concentration parameter, formation redshift, spin parame- 
ter and shape and found that halo mass is the only property 
that correlates significantly with local environment. It is im- 
portant to note that the variation of the halo mass function 
in different environments, ie., the fact that high mass haloes 
are under-represented and over-represented in low and high 
density regions, respectively, suggests that any apparent de- 
pendence of halo properties on the environment, could be a 
consequence of the dependence of these properties on halo 
mass. Therefore, one needs to disentangle the two dependen- 
cies and to perform any environmental dependence study as 
a function of halo mass as well. 

A large number of recent studies on the environmen- 
tal effects on a variety of halo properties, like halo shapes, 
spin, alignments velocity dispersion, and for different mass 
halo ranges, have been presented (eg. Faltenbacher et al. 
2002; Einasto et al. 2003; Ragone et al. 2004; Einasto et 
al. 2005; Avila-Reese et al. 2005; Hopkins, Bahcall & Bode 
2005; Basilakos et al. 2006; Plionis et al. 2006; Altay et 
al. 2006; Maulbetsch et al. 2006; Aragon-Calvo et al. 2006; 
Hahn et al. 2006). However, results of different studies are 
not always in agreement with each other, a fact that could 
be due to different quantifications of the environment or due 
to different analysis tools. For example, Einasto et al. (2005) 
found that group and cluster size DM haloes, in high den- 
sity regions, have smaller eccentricities (are more spherical) 
than in low density regions, while Kasun & Evrard (2005) 
have found no such dependence. Similarly, Avila-Reese et al. 
(2005) and Hahn et al (2006) have found a dependence of 
galaxy-size DM haloes shapes on environment but again no 
such obvious dependence for larger DM haloes. 

An interesting property to study as a function of en- 
vironment, is the internal alignment between the principal 
axes of the shape and velocity anisotropy ellipsoids, which 
can be considered as an indication of relaxation in a system 



where the shape is supported by internal velocities (eg. Tor- 
men 1997). Kasun & Evrard (2005) and Allgood et al. (2006) 
found for cluster-size DM haloes a good such alignment, al- 
though no investigation in different environments has been 
reported. Furthermore, the external alignment between DM 
halo axes or angular momentum and the orientation of the 
filament in which the halo is embedded is of interest. Bailin 
& Steinmetz (2005) found a very strong tendency for the 
halo minor axis to lie perpendicular to the large scale fila- 
ment, but a much weaker tendency for the major axis to be 
oriented parallel to it. They also found that the group and 
cluster size halo angular momenta lie perpendicular to the 
large scale filaments while that of galaxy-size haloes tend 
to lie parallel to them. This suggests that group-size DM 
haloes acquire most of their angular momenta from mergers 
along the filament direction. Avila-Reese et al. (2005), in 
turn find a decreasing alignment signal between minor axis 
and angular momentum of galaxy-size DM haloes going from 
overdense to underdense regions. 

In this work we attempt to give a new insight in the 
behavior of group-size DM halo properties (shape, velocity 
dispersion, internal and external alignments) as a function of 
environment, taking special care to disentangle their corre- 
lation with halo mass, as mentioned before. We also divide 
our sample of group-size haloes according to the amount 
of substructure that they have, in order to infer if merg- 
ers and/or gravitational tidal interactions play a significant 
role in shaping the DM halo morphological and dynamical 
properties. 

The outline of this paper is as follows. In section 2 we 
describe the numerical simulation method, the halo iden- 
tification procedure and the research methodology that we 
will follow. In section 3, we present the methods for the 
computation of the shape and velocity tensors, the angular 
momentum and alignment measures and finally we present 
a thorough study for the quantification that we use to deter- 
mine the halo substructure. In section 4 we study the depen- 
dence of halo shape and dynamics on environment while in 
section 5 we present the corresponding study of internal and 
external halo alignments. Finally, we summarize our results 
and draw our conclusions in section 6. 



2 NUMERICAL DATA & RESEARCH 
METHODOLOGY 

The numerical simulations used in this work were performed 
using the GADGET2 code (Springel 2005) with dark matter 
only. This parallel code was ran in a Beowulf cluster with 
32 Intel Xeon processors (3.06 GHz). The cosmological pa- 
rameters used correspond to a flat cosmological model with 
a non- vanishing cosmological constant (ACDM): Q. m = 0.3, 
£Ia = 0.7, as = 0.9, h = 0.72, where Q, m and Qa are the 
present day matter and vacuum energy densities in units 
of the critical density, as is the present linear rms ampli- 
tude of mass fluctuation in spheres of 8ft" 1 Mpc and h is 
the Hubble parameter in units of 100km s _1 Mpc _1 . The 
initial conditions are generated with the GRAFIC2 pack- 
age (Bertschinger 2001), which also computes the transfer 
function as described in Ma & Bertschinger (1995). 

The main simulation was run in a cube of size L — 
500/i _1 Mpc, using 512 3 particles. The particle mass is ~ 
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7.7 x lO lo /i -1 M0 and the force softening length is e = 
100/i _1 kpc. Individual particle time-steps are chosen to be 
proportional to the square r oot of t he softening length over 
the acceleration a: AU = y/2rje/\a\. We set the dimension- 
less parameter which controls the accuracy of the time-step 
to be 7] = 0.02. 

The haloes were identified using a FOF algorithm with 
a linking length I — 0.17 times the mean inter-particle sepa- 
ration. Given the proposes of this work, we only use haloes 
with at least 130 particles, ie., with masses greater than 
1O 13 /i _1 M0. Note that this halo finder does not identify (for 
a given linking length) sub-haloes belonging to larger par- 
ent haloes. For the purpose of our study, however, we will 
consider haloes with and without substructure (see further 
below), which in effect correspond to those haloes with or 
without relatively massive sub-haloes. 

The resulting sample of ~ 58000 haloes was split in 
two subsamples: haloes with masses M > 2 x 1O 14 /i~ 1 M0 
are considered as clusters (in total 1598 haloes), whereas 
haloes in the range 10 13 fr _1 M Q < M < 2 x 10 14 h _1 M o are 
considered as groups (56699 haloes). 

In order to investigate the role that environment plays 
in determining halo properties, we find for each halo the 
distance to its nearest cluster and divide the halo sample in 
three subsamples according to this distance (r c i ustcr ): 

• Small distance subsample: r c i ustcr < 7/t _1 Mpc (H0-7 
hereafter, ~ 8% of the group sample). 

• Intermediate distance subsam- 
ple: 10/i _1 Mpc < r c i ustor < 17/i _1 Mpc (H10-17 hereafter, 
~ 21% of the group sample). 

• Large distance subsample: 30/i _1 Mpc < r c i ustC r < 
50/i -1 Mpc (H30-50 hereafter, ~22% of the group sample). 

The left panel of Figure 1 shows the fraction of haloes in each 
one of the previously defined subsamples. In order to ensure 
that r c iu S tor is defining accurately the environment, we also 
compute for our haloes at different distances from the clus- 
ters the corresponding density contrast 5(r) = p(r)/p — 1, 
where p(r) is the density in a sphere of radius 8/i _1 Mpc 
around the halo centers and p is the mean matter back- 
ground density. Results are shown in the right panel of Fig- 
ure 1 from which it is obvious that indeed the halo distance 
to its nearest cluster is related to the overdensity in which 
the halo is embedded. 

Given the significant effect that mergers and interac- 
tions can have on the shapes and alignments of haloes, and 
the fact that in overdense regions the halo mass function 
is skewed toward the high-mass end (Lemson & Kauffmann 
1999), we will present our results as a function of their: 

(a) halo mass, 

(b) environment, determined by r c i ustcr , and 

(c) halo dynamical state, determined by the Dressier & 
Shectman (1998) method (see next section). 

Therefore, any dependence of the halo properties on 
environment will be disentangled from the mass function 
effects, and will not be attributed to the overabundance of 
high mass haloes in overdense regions. 

For the purpose of testing the robustness of our sub- 
structure determination procedure to variations of the sim- 
ulation resolution and box size, we also run (a) one sim- 
ulation with the same resolution as the main simulation 
but in a 8 times smaller box, ie., evolving 256 3 particles 
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Figure 1. Left panel: Fraction of haloes in the H0-7 (vertical 
dashed histogram), H10-17 (empty histogram) and H3o_5o (hor- 
izontal dashed histogram) sample. Right panel: The median den- 
sity contrast 5(r) computed in spheres of radius 8h~ 1 M.pc, as a 
function of r c i ustcr distance. Error bars represent the 33% and 
67% quantiles of the corresponding distribution. 

in a L = 250/i~ 1 Mpc side box (LR hereafter), and (b) a 
higher resolution simulation obtained by re-simulating the 
L — 125/i -1 Mpc central box of the former with 256 3 par- 
ticles and e = 50/i _1 kpc (HR hereafter), reaching a parti- 
cle mass resolution of 9.7 x 1O 9 /i _1 M0. We identify haloes 
in both the LR and HR simulations using a linking length 
I = 0.17 times the mean inter-particle separation as in our 
main simulation. 

3 DETERMINATION OF HALO SHAPE, 
ALIGNMENTS & DYNAMICAL STATE 

3.1 Parameter Definition 

The shape of haloes, modeled as ellipsoids, is determined by 
diagonalizing their inertia tensor: 

I%j ~ ^ ^ %i,n%j,n , (1) 
JV 

where N is the number of particles in the halo and Xi, n is 
the i th component of the position vector of the n th particle 
relative to the halo center. The principal axes of the fitted 
ellipsoid (a, b, c with a > b > c) are related to the square root 
of the eigenvalues of the inertia tensor. The corresponding 
eigenvectors provide the directions of the principal axes of 
the fitted ellipsoid. 

Similarly, velocity moments are obtained by diagonaliz- 
ing the velocity anisotropy tensor: 

Vij = y]vi, n Vj,n , (2) 
N 

where v i>n is the i th component of the velocity vector of 
the n th particle relative to the halo center of mass velocity. 
Note that a vc i > fe vc i > c vc i will denote the major, middle 
and minor axes of the velocity ellipsoids, respectively. 

We compute the specific angular momentum of each 
halo containing N particles as: 

L =^E r,XVl ' ( 3 ) 

N 

where ri and Vi are the position and velocity vectors of the 
particle i relative to the halo center of mass. 
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The various alignments between different pairs of vec- 
tors, representing either the principal axes of the halo den- 
sity and velocity ellipsoid, the halo angular momentum or 
the direction to a neighboring cluster halo, will be estimated 
by the mean of the distribution of | cos(#)|, where 9 is the 
angle between the directions of any two vectors, vi and v 2 , 
we are interested in. Therefore, 



COS(0) = Vl • V2 



(4) 



Perfect alignment and ant i- alignment correspond to 
|cos(0)| = 1 and 0, respectively, whereas for the random 
three-dimensional case the expected distribution mean value 
is (|cos(0)|) = 0.5. 

Finally, we use the Dressier & Shectman (1998) algo- 
rithm to estimate the amount of substructure in haloes. 
Briefly, this method determines the mean local velocity 
(vioc) and the local velocity dispersion ui oc of the nearest 
n neighbors from each halo particle i and compares them 
with the mean velocity, (V), and the velocity dispersion, 
a, of the whole halo of N particles, defining the following 



^ 2 = -[({v loc )-{V)) 2 + (a loc -a) 2 ] 

a 



where 



and 



En( V l°c - Vi) 

n-\ 



Ejv((V)-Vi 
N — 1 



(5) 



(6) 



(7) 



A quantification of the substructure present in a halo is 
given by the so-called A- deviation, which is the sum of the 
individual Si's over all halo particles N: 



A = 



T,N S i 



N 



(8) 



The larger the A-deviation the stronger is the halo sub- 
structure. This statistic depends on the number of nearest 
neighbors n which is used in the analysis, and as we verified 
on the number of particles used to resolve a halo as well. 

We have computed the A-deviation using two different 
values of n: (a) n = 25 as in Knebe & Miiller (1999) and (b) 
n = N 1 / 2 as in Pinkney et al (1996) and find similar results. 



3.2 Random and Systematic Parameter 
Uncertainties 

We investigate the uncertainty introduced by resolution ef- 
fects in the determination of our morphological and dynam- 
ical halo parameters. The fact that low mass haloes are 
resolved by a smaller number of particles with respect to 
higher mass haloes would inevitably create a random or pos- 
sibly even a systematic deviation from their nominal values. 
To investigate these uncertainties we use a procedure similar 
to that of Avila-Reese et al. (2005) . 

We perform 100 realizations of each massive halo (hav- 
ing more than 5000 particles) which we resolve selecting 
randomly the same number of particles as that of the low- 
est mass haloes used in our analysis (ie., 130 random parti- 
cles). For each realization we then compute the halo c/a axis 
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Figure 2. la error distribution of the various dynamical and 
morphological parameters of haloes with more than 5000 particles 
but sampled with only 130 random particles (100 realizations are 
used). The distributions of the halo axial ratios, velocity disper- 
sions, velocity-shape major axes mis-alignment angles and minor 
axis-angular momentum mis-alignment angles are shown in the 
Top left, Top right, Bottom left and Bottom right panels, respec- 
tively. 



ratio, velocity dispersion, velocity-shape major axes mis- 
alignment angle, a-a ve i, and minor axis-angular momentum 
mis-alignment angle, c • L. 

The distributions of the la deviations from their nom- 
inal value, plotted in Figure 2, have mean values of ~ 
9%, 3.5%, 16% and 30%, respectively for the c/a axis ra- 
tio (top left panel), velocity dispersion (top right panel), 
velocity-shape major axes mis-alignment angle (bottom left 
panel) and minor axis-angular momentum mis-alignment 
angle (bottom right panel). It is evident that the uncertain- 
ties are quite small, especially of the halo velocity dispersion, 
except for the mis- alignment angle, fi • L. 

It is also possible that resolution effects do not only 
introduce a random error on the nominally defined shape 
and dynamical halo parameter. For example, for the case of 
the substructure index, A, a pronounced trend is apparent 
with (A) increasing with halo mass. In Figure 3 (left panel) 
we show results based on the n = N 1 ^ 2 case (filled circles & 
solid line). There is an apparent monotonic increase of (A) 
with halo mass. Furthermore, in the right panel of Figure 3 
we show with the solid line, the percentage of haloes with 
A-deviation higher than the mean of the A distribution of 
all haloes ((A) = 0.98). Again, we see a monotonic increase 
of the fraction of haloes having substructure as a function 
of mass, with the most massive haloes appearing all to be 
substructured. These results create a suspicion that they 
could be due to the lower resolution with which the low 
mass haloes are resolved. To investigate the resolution issue 
we perform two tests: 
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Figure 3. Left panel: The dependence of the mean A-deviation 
substructure index on halo mass. The filled circles, connected 
with the solid line, represent the mean A-deviation computed 
using all the halo member particles, whereas the corresponding 
values choosing to have the same halo resolution, ie., ~ 130 ran- 
dom particles per halo (A ran ), is shown as open circles connected 
with the dashed line. The estimated bootstrap uncertainties are 
smaller than the size of the symbols. Right panel: The percentage 
of haloes with A-deviation higher than the mean value (ie., haloes 
with substructure), as a function of halo mass. The solid line and 
filled symbols correspond to haloes selected using the global dis- 
tribution of A-deviations, estimated using all halo particles, irre- 
spective of the halo mass (variable halo particle resolution). The 
dashed line and open symbols correspond to haloes selected using 
the global distribution of A ran -deviations, estimated using ~ 130 
particles per halo, irrespective of the halo mass (same halo par- 
ticle resolution). While the dotted line shows the percentages of 
haloes with A > (Aj), evaluated within the i th bin of halo mass 
and using all halo particles. 



• we compare the A-deviation index for the matching 
haloes of our HR and LR simulations, and 

• we recompute A for all mass haloes, but using only 130 
randomly selected particles per halo (i.e. the same number 
resolution as in the smaller haloes). 

Regarding the first test we select those pairs of haloes 
which match, in position and mass, in both the LR and 
HR halo samples. The masses of these matching haloes are 
allowed to differ only by 5%, so as to ensure that we will 
compare properties of the same haloes with only difference 
their resolution (~ 8 times more particles in the HR match- 
ing haloes). We divide the matching haloes sample in three 
subsamples according to their masses, which are compared 
in the left panel of Figure 4. In the right panel we com- 
pare their corresponding A-deviation values and in all cases 
we find that A is significantly larger when computed in the 
higher resolution haloes, a fact which is further enhanced for 
the more massive haloes. This result verifies our suspicion 
that resolution effects could be the cause of the monotonic 
increase of A with halo mass. 

We now continue with our second test and derive for 
each halo a new A-deviation index (A ran ) computed by us- 
ing the same number of particles (130), randomly selected, 
in each halo independent of its mass. In this way we impose 
the same particle resolution on all haloes. Note that we use 
as A ran the average over many realizations of the random 
particle selection process. The results of this procedure show 
that A ran is systematically smaller than when using all the 
particles in the haloes, and increases very weakly with halo 




Figure 4. Left panel: Comparison of the masses of matching LR 
and HR haloes: open circles, triangles and squares correspond to 
haloes with masses: 10 13 < M < 3xlO 13 fe _1 M , 3xl0 13 < M < 
6 x 10 13 /i _1 Mq and 6 x 10 13 < M < 1O 15 /i _1 M , respectively. 
Right panel: The A-deviation comparison of the matched LR and 
HR haloes. 

mass, as shown by the dashed line in both the left and right 
panels of Figure 3. Again these results indicate the impor- 
tance of resolution effects in quantifying the amount of halo 
substructure. 

From our previous study we have realized that the whole 
distribution of A-deviations shifts to higher values as a func- 
tion of resolution. This translates to a shift at higher values 
of the A-deviation distribution as a function of halo mass, 
within the same simulation. 

In the right panel of Figure 3 we also show, as the dot- 
ted line, the percentage of haloes as a function of halo mass, 
having their A estimated using all halo particles, but then 
selecting those with A larger than the mean, (Abin), of the 
distribution within each bin of halo mass. It is evident from 
this plot that this case and the case based on A ran give 
equivalent results. Had we not taken into account the res- 
olution effects we would have erroneously concluded that 
almost all massive haloes have strong substructure and that 
the opposite was true at the low halo mass end. 

Given the computation of the substructure A-deviation 
index is more robust if all halo member particles are consid- 
ered, we will refer from now on to haloes having substruc- 
ture as those with A > (Abin) in the specific mass range bin 
which they belong. 

The fraction of haloes with substructure, as defined be- 
fore, are ~ 45% ± 5% for the H _ 7 haloes and 35%±3% for 
the H10-17 and H30-50 haloes, with only a weak dependence 
on halo mass. If, however, we select haloes nearer to mas- 
sive clusters (ie., < r cluster < 4 h" 1 Mpc), the fraction of 
low-mass haloes (M < 2 x 10 13 H^ 1 Mq) with substructure 
grows to 65% ±3%. 

Note that haloes that went through a recent merger 
will have a higher A-deviation value with respect to those 
that are either isolated or had no recent merger event, thus 
having more time to virialize. 

In Figure 5 we present the dependence of the substruc- 
ture index on the environment and halo mass. Solid and dot- 
ted lines stand for the ratios of A in the H0-7 and H10-17 
samples, normalized to the most distant sample (H30-50). As 
expected, we find haloes in the vicinity of massive clusters 
(solid line) to have larger A values (for their mass range) 
with respect to distant haloes (dotted line), presumably due 
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Figure 5. A-deviation ratios as a function of halo mass and envi- 
ronment. The solid line represents the ratio: (Ah _ 7 )/{Ah* 30 _ 50 ) 
and the dotted line the ratio: (Ajj 10 _ 17 )/(Ah 30 _ 50 )■ Error bars 
are based on the propagation of the individual A-deviation un- 
certainties. 



to the higher merging rate and due to the stronger tidal 
field, found around overdense regions. 



4 ENVIRONMENTAL EFFECTS ON HALO 
SHAPES & DYNAMICS 

4.1 Halo Shape - Mass correlation 

In the ACDM cosmology, the dependence of shapes on DM 
halo mass has been well established in many recent studies 
with more massive haloes being less spherical, ie., having a 
lower axis ratio, c/a (eg. Bullock 2002; Jing & Suto 2002; 
Kasun & Evrard 2005; Gottlober & Turchaninov 2006; All- 
good et al. 2006; Paz et al. 2006; Maccio et al. 2006; Bett 
et al. 2006). This can be explained considering that in the 
hierarchical clustering of CDM haloes, smaller mass haloes 
form earlier on average than massive ones, and thus they 
have more time to evolve, virialize and become more spheri- 
cal. We should also note that including baryonic physics has 
a significant effect on the shapes of haloes (eg., Kazantzidis 
et al. 2004). In this section we investigate whether this 
trend changes when considering groups in different envi- 
ronments. Such a difference has been noted by Avila-Reese 
et al. (2005) between galaxy size haloes found in clusters 
and in voids, and by Hahn et al. (2006) between haloes 
found in clusters and in filaments but only for small halo 
masses (M < 2 x 10 12 M Q ). We do not probe this mass 
range and therefore our analysis concentrates only on larger 
mass haloes, typical of groups and poor clusters of galaxies. 

Figure 6 (left panel) presents our results for the three 
r cluster subsamples (H -7, H10-17 and H30-50), with error 
bars computed using the bootstrap re-sampling technique. 
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Figure 6. Left panel: Mean axial ratios {c/a) as a function of halo 
mass and environment. Dotted, short dashed and long dashed 
lines correspond to halo samples of different cluster-group dis- 
tances (Ho— 7, H10— 17 and H30— 50, respectively). The standard 
deviation of the axial ratio distribution for the different halo sub- 
samples and halo masses is between 0.11 and 0.12. Error bars were 
calculated using the bootstrap re-sampling technique. The solid 
line denotes Allgood et al. (2005) fit. Right panel: Mean axial ratio 
(c/a), as a function of halo mass, only for the Ho_7 haloes (dot- 
ted line). Squares and triangles correspond to Ho -7 haloes with 
A-deviation values lower and higher than (Abm)i respectively. 



In all cases the trend we find is in accordance with the well 
known mass-shape relation, with (c/a) increasing with de- 
creasing halo mass. However, it is shifted toward more spher- 
ical axial ratios when considering haloes nearer to massive 
clusters. This can be explained by the fact that haloes in 
high-density environments are formed earlier than haloes, of 
the same mass range, in low-density environments (eg. Sheth 
& Tormen 2004) , giving the former more time to evolve, re- 
lax and hence become more spherical (Avila-Reese et al. 
2005). This is in agreement with haloes at higher redshifts 
being more elongated than present day equal mass haloes 
(eg. Allgood et al. 2005). 

Another representation of our results is shown in Fig- 
ure 7 where we plot the c/a frequency distribution of well re- 
solved haloes (ie., those with 3.8 x 10 13 < M < 10 14 /i _1 M a ). 
The means of these distributions are 0.50, 0.48 and 0.47 
for the Ho-7(dashed histogram), H10-17 (dotted histogram) 
and H30-50 (empty histogram) samples, respectively, while 
a Kolmogorov-Smirnov two-sample test shows them to be 
different at a very high significant level. 

We now consider only the H0-7 haloes and compute the 
shape-mass relation but separating haloes with high and low 
A-deviation (substructure). Results are plotted in the right 
panel of Figure 6, where squares correspond to haloes with 
A < (Abin) and triangles to haloes with A > (Abin). The 
shape-mass relation is basically maintained in both subsam- 
ples although: (a) it is much shallower for haloes in the low 
A-deviation subsample, and (b) it is shifted toward lower 
(c/a) for haloes in the high A-deviation subsample. This 
result is indeed expected if to consider that more virial- 
ized systems tend to be more spherical. Those system with 
high level of substructure, which are dynamically younger 
systems, have interrupted their virialization process due to 
some recent merger event and therefore have a more elon- 
gated shape than systems with no sign of substructure. This 
latter behavior is present in the whole range of considered 
masses and also in the H10-17 and H30-50 subsamples. 
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Figure 7. Axial ratio (c/a) distributions for haloes in the mass 
range 3.8 X 10 13 - 10 14 /i~ 1 M Q . Shaded, dotted and empty his- 
tograms correspond to the Ho— 7, H10— 17 and H30-50 halo sam- 
ples, whose means are 0.50, 0.48 and 0.47, respectively. The la 
dispersion is in all cases ~ 0.1 

In order to discard the possibility that discreteness ef- 
fects could impose the (c/a)-Mass correlation, we recom- 
puted the halo shapes but using realizations of only 130 
randomly selected particles per halo, so as to resemble the 
resolution of the less massive haloes (see also Paz et al. 2006 
for discreteness related effects) . We indeed recover the same 
(c/a)-Mass trend and thus we verify that it is not imposed by 
the variable resolution with which the different mass haloes 
are resolved. However, there is a shift toward less spheri- 
cal values when using the common halo resolution of 130 
randomly selected particles per halo. The apparent curva- 
ture, toward lower c/a values and at the low-mass end of 
the (c/a)-Mass relation, apparently disappears when using 
the common halo resolution and thus it should probably be 
attributed to the variable halo resolution. 

4.2 Halo Velocity Dispersion-Mass correlation 

In Figure 8 we present the halo velocity dispersion-mass cor- 
relation. Velocity dispersions were computed using formula 
7. Such a correlation is expected from to the virial theorem. 
In order to investigate the possible influence of the environ- 
ment on this relation and hence on the reliability of using 
the virial theorem to estimate halo masses, we present in the 
left panel of Fig. 8 results for the three halo r c i ustcr (H0-7, 
H10-17 and H30-50) subsamples and for an extra subsample 
with r c i us tcr < 4/i _1 Mpc (dot-dashed line). 

As expected from the virial relation we find that larger 
mass haloes have higher velocity dispersions, for all the con- 
sidered subsamples. However, there is a shift toward higher 
velocity dispersions of low-mass haloes found near clusters 
with respect to those found further away. This trend is 
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Figure 8. Mean velocity dispersion (a) as a function of halo 
mass. Left panel: Results for the usual three r c i ustcr subsam- 
ples (H0-7, H10— 17, H30-50), with line types already defined in 
Figure 6, while the dot dashed line corresponds to haloes with 
r cluster < 4 /i -1 Mpc (the H10— 17 and H30— 50 results are iden- 
tical). Bootstrap errors are small, typically ~ 10 — 20 km/scc. 
Right panel: Results only for the Ho— 7 halo subsample Squares 
correspond to haloes with A < (A^in) and triangles to haloes 
with A > (A bin ) . 

stronger the nearer the low-mass halo is found to the cluster. 
However, for halo-cluster distances > 10/i _1 Mpc there is no 
effect whatsoever. 

To investigate whether the halo dynamical state relates 
to the halo velocity dispersion-mass correlation, we divide 
the H0-7 subsample, as in the previous section, to those 
with and without substructure. Results are shown in the 
right panel of Figure 8, where the dotted line stands for 
all the H0-7 haloes, squares and triangles for haloes with 
A < (Abin) and A > (Abi n ) , respectively. There is indeed 
a dependence of the mass- velocity dispersion correlation on 
the amount of halo substructure but only for haloes with 
masses < 5 x 10 13 hT 1 Mq, which are found to have a 
larger mean velocity dispersion than haloes with no or low- 
levels of substructure. Note that Evrard et al. (2007) find 
similar results for what they call satellite haloes. Moreover, 
we find that haloes with a low A-deviation index behave 
similarly as haloes in the H10-17 and H30-50 samples (seen 
in the left panel of Figure 8), which show no dependence 
of the mean velocity dispersion on the presence or not of 
halo substructure. Note, however, that a slight shift toward 
higher velocity dispersion values is present also in high mass 
haloes but only if choosing those haloes with extremely high 
A-deviation index. 

4.3 Partial Conclusions 

These results, concerning the dependence of halo shapes and 
velocity dispersion on the halo dynamical state, give new in- 
sights in our understanding of halo formation and evolution. 

Although the general expectation is that low-mass 
haloes in high density environments formed earlier and thus 
should be relatively more virialized with respect to similar 
mass haloes in low-density regions, we have found a rela- 
tively high fraction of dynamically young and active haloes 
near massive clusters (;> 45%). These haloes have in gen- 
eral a higher velocity dispersion (more evident at the low 
mass end) and a lower {c/a) ratio with respect to similar 
mass virialized haloes. The high level of substructure of these 
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haloes is probably because they are continuing to grow via 
mergers in the anisotropic outskirts of massive haloes (eg. 
West 1994; Maulbetsch et al. 2006), although their dynam- 
ical state could also be affected from the strong tidal field 
imposed by their local high density surrounding, while both 
cases imply a lower halo sphericity, as observed. The mergers 
as the most possible cause for the increase of the halo veloc- 
ity dispersion is in agreement with Faltenbacher et al. (2006) 
who find for an equal mass merging event (progenitors with 
masses ~ 1 x 10 14 h^ 1 MQ) an oscillatory behaviour of the 
velocity dispersion (among other properties). After the re- 
laxation of the new system, the velocity dispersion is slightly 
larger, but it changes substantially during the event. 

Now, higher mass haloes with a high A-deviation (sub- 
structure) in the vicinity of massive clusters also appear to 
be of lower sphericity although, and contrary to the low- 
mass halo case, their velocity dispersion does not show any 
significant deviation from that of the more virialized high- 
mass haloes. This could be explained if typically the merger 
events, which alter the higher mass halo shape, are due to 
relatively lower mass haloes which although affect the overall 
shape, they affect less the dynamical structure of the high 
mass halo, which is dominated by the main gravitational 
potential of the high-mass halo itself. It could also imply a 
faster "re-accommodation" of the velocity field with respect 
to the density field, in the relatively deep principal halo po- 
tential well. 



5 ENVIRONMENTAL EFFECTS ON HALO 
ALIGNMENTS 

5.1 Internal Halo Alignments 

In section 1 we mentioned that there is a significant signal 
of alignment between the principal axes of the shape and 
the velocity anisotropy tensors, indicating that most of the 
haloes have their shapes supported by the velocities of their 
member particles (eg. Tormen 1997; Kasun & Evrard 2005) . 

Here we investigate whether there is any environmen- 
tal dependence on such an alignment effect. To this end we 
compute the mean absolute cosine of the angle between the 
major axes of the mentioned tensors, (|a • a vc i|), as a func- 
tion of the halo mass for the three halo subsamples (H0-7, 
H10-17 and H3o-so). 

In the left panel of Figure 9 we show the (|a • a ve i|)- 
mass correlation for the H0-7 (dotted line), Hio_i7 (dashed 
line) and H30-50 (long dashed line) samples. In all cases 
there is a good signal of alignment between the shape and 
velocity ellipsoid principal axes, specially for haloes at large 
distances from massive clusters (H10-17 and H30-50 subsam- 
ples), while within each subsample the alignment is stronger 
for the higher mass haloes. We have verified that this is not 
due to the variable resolution with which the different mass 
haloes are sampled. 

Haloes at large distances from massive clusters have 
their velocity and shape better correlated probably because 
they are less tidally disrupted than in the high density en- 
vironment of the cluster, either by the cluster itself or/and 
by the local overabundance of lower mass haloes found in 
such environment. Even more so in the high-mass halo end 
probably because interactions and merging with lower mass 
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Figure 9. Left panel: the direction cosine of the major axes of the 
shape and velocity ellipsoids ([a-a ve i|) as a function of halo mass 
for the Ho_7, Hio_i7 and H3o_so halo subsamples. Line styles 
are as in Figure 6. Error bars were calculated using the bootstrap 
re-sampling technique. Right panel: the corresponding (|a-a ve i|)- 
mass correlation only for the Ho- 7 sample, split in those haloes 
with high level of substructure, A > (A bin ) (triangles), and those 
without substructure, A < (A^in) (squares). 

haloes can disturb minimally the phase-space of these higher 
mass haloes. However, near the massive cluster (H0-7 sub- 
sample), the stronger halo-cluster gravitational interactions 
affect significantly the halo phase-space and for this reason 
we observe a general decrease of the value of (|a • a ve i|) for 
all halo masses. 

This interpretation could be supported if those haloes 
with a high level of substructure showed even less aligned 
orientations. Indeed this is the case, as can be seen in the 
right panel of Figure 9, where we present for the H0-7 
subsample the (|a • a ve i|)-mass correlation but split be- 
tween haloes with high level of substructure, A > (A b m) 
(triangles), and haloes with no substructure, A < (Abin) 
(squares). The former haloes show a weaker alignment, as 
anticipated, suggesting that strong interactions and mergers 
introduce scatter in the phase space of these systems. 

The segregation between haloes with and without sub- 
structure, seen in the right panel of Figure 9, is also present 
in the more distant halo samples (H10-17 and H3o-so). The 
same interpretation, given before, of the difference between 
equal mass haloes with and without substructure, holds for 
these samples as well. Furthermore, the slightly better align- 
ment seen for the more massive haloes is due to their deeper 
potential wells, which inevitably creates a better alignment. 

Another internal alignment effect that we address is 
that between the directions of the angular momentum vec- 
tor and the minor axis of the mass distribution, c • L|. It 
has been found that the angular momentum is most often 
aligned with the minor axis and perpendicular to the ma- 
jor axis (eg. Dubinski 1992; Bailin & Steinmetz 2005; Bett 
et al. 2006). An environmental dependence has also been 
found by Avila-Rcese et al. (2005) for galaxy-sized haloes, 
with a higher alignment signal in underdense regions (see 
also Hahn et al. 2006 for angular momentum orientations 
with large-scale structures). 

Our results are shown in Figure 10, where we also find 
such an alignment signal, although it is obvious that due to 
noise we are unable to detect different trends in the three 
halo subsamples. Note also that the amplitude of our align- 
ment signal is significantly less than that found by Bailin 
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Figure 10. Left panel: the direction cosine between the angular 
momentum and the minor axis ( | c ■ L | ) as a function of halo mass 
for the Ho_7, Hio_i7 and H3o_so halo subsamplcs. Right panel: 
the corresponding (|c ■ L|)-mass correlation only for the H0-7 
sample, split in those haloes with high level of substructure, A > 
(Abi n } (triangles), and those without substructure, A < (At,; n ) 
(squares). 
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Figure 11. Left panel: Correlation between halo-cluster direction 
and halo major axis orientation ((|dg n ■ a|)) as a function of halo 
mass. Long dashed, short dashed and dotted lines correspond to 
different bins of cluster-group distances as in Figure 6. Error bars 
were calculated using the bootstrap re-sampling technique. Right 
panel: Similar but only for the H0-7 sample, split into haloes with 
and without significant substructure, A > (A bin ) (triangles) and 
A < (Abi n ) (squares), respectively. 



& Steinmetz (2005), most probably because these authors 
define shapes using the reduced moment of inertia tensor, 
which weights strongly the inner parts of haloes, as well as 
because they choose to analyse only haloes of which both 
their small axis and angular momentum orientation have 
small uncertainties. Furthermore, one should keep in mind 
that, due to resolution effects (see section 3.2), the intrin- 
sic uncertainty of this alignment measure is quite large for 
low-mass haloes. 

Returning to our results we do find a systematic, al- 
though weak, trend of a better alignment for the H0-7 sub- 
sample (left panel of Figure 10), which is in the opposite 
direction than the results of Avila-Reese et al. (2005) based 
on galaxy-sized haloes. However, our haloes are much larger 
(groups and poor cluster size) and this could well be the 
reason of the apparent discrepancy. Furthermore, we are in 
general agreement with the recent results of Aragon-Calvo 
et al. (2006). 

Moreover, haloes with high level of substructure seem 
to have L and c better correlated (right panel Figure 10), a 
fact which is true for all halo subsamples (H0-7, H10-17 and 
H3o-5o)- This should be partly attributed to the fact that 
haloes with substructure are more elongated than relaxed 
haloes and thus they have both their angular momentum 
and minor axis vectors better defined. Furthermore, this re- 
sult also implies that mergers probably affect significantly 
the angular momentum of haloes (eg. Vitvitska et al. 2002), 
which gain part of their angular momentum from mergers 
preferentially occurring along the plane defined by the major 
and median axes. 



5.2 External Alignments 

It has been shown that the orientation of the halo major axis 
is strongly correlated with the direction from which the last 
major merger event occurred (van Haarlem & van de Wey- 
gaert 1993). Therefore, it should be expected to find a corre- 
lation between halo major axis orientation and the direction 
denned by the halo-cluster distance, which in turn should in- 
dicate the orientation of the filament. Such alignment effects, 



among relatively massive haloes, have been found to be par- 
ticularly strong and extending up to > 100 h' 1 Mpc (eg., 
Faltenbacher et al. 2002; Kasun & Evrard 2005; Hopkins, 
Bachall & Bode 2005). 

Here we also wish to calculate the alignment between 
the direction of each halo to their nearest cluster (d gn ) 
and the halo major or minor axis orientation ((|d gn • a|) 
and (|d g h • c|), respectively) for all three (H0-7, H10-17 and 
H30-50) halo subsamples. The left panel of Figure 11 shows 
the case for the major axis alignment while the left panel 
of Figure 12 shows the corresponding minor axis alignment 
case. As expected, the former alignment effect is strong, 
more so for haloes found near their clusters and for the high 
mass haloes. As in Bailin & Steinmetz (2005) we find that 
the halo minor axes are in general anti-aligned (perpendicu- 
lar) to the filament direction, again more so for haloes found 
near their clusters and for high mass haloes. Now dividing 
our Ho -7 subsample into those haloes with and without sub- 
structure (right panels of Figures 11 and 12), we find very 
interesting results: 

• relatively virialized haloes, having no significant sub- 
structure, show a strong tendency for major axis alignment 
(and minor axis anti-alignment) with the direction to their 
nearest massive cluster implying that they retain strong 
memory of the initial anisotropic distribution from which 
they accreted matter (eg. van Haarlem & van de Weygaert 
1993), 

• haloes with high level of substructure show a similar, 
although relatively weaker, alignment effect, which appears 
to be in disagreement with what one would naively expect 
given that mergers happen preferentially along the filaments. 
However, once a merger has happened, non-linear gravita- 
tional effects take place and until the merged structure re- 
laxes, it may well appear less aligned with the filament ori- 
entation. Specially, low mass haloes, being also small in size, 
interacting with other neighboring haloes in the high den- 
sity surroundings of a massive cluster, could be relatively 
more affected by local gravitational effects which may not 
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Figure 12. Left panel: Correlation between cluster-group direc- 
tion and halo minor axis orientation (ddgh ■ 6|)) as a function of 
halo mass. Right panel: Similar but only for the H0-7 sample, 
split into haloes with and without significant substructure. 



necessarily reflect the large-scale anisotropic distribution of 
matter in the filament. 

An interesting complication in the above interpretation 
is that the halo angular momentum and minor axis is better 
aligned in high A-deviation haloes (see section 5.1) and if di- 
rectional mergers are responsible for such an alignment, be- 
ing more frequent in the filament direction, then one might 
have expected their major axes to be more aligned with the 
halo-cluster direction than for virialized haloes (small A- 
deviations). However, what the angular momentum- minor 
axis halo alignment actually implies, in the above picture, 
is that the direction of the merger is in the plane defined by 
the major and medium axes and not necessarily along the 
major axis, which is not in contradiction with the above. 



6 DISCUSSION & CONCLUSIONS 

We have used haloes identified with a FOF algorithm in 
a dark matter only ACDM simulation to study the de- 
pendence of the shape, dynamical state and various align- 
ments of group size haloes on the halo environment. The 
smallest haloes analysed have at least 130 particles (M ;> 
10 13 /i _1 Mq). In order to investigate if there is some influ- 
ence of the environment on the properties of haloes, we split 
the group sample in three subsamples according to their dis- 
tances to the nearest massive cluster halo. We have also in- 
vestigated the Dressier & Shectman (1996) algorithm, used 
to determine whether a halo has a high level of substruc- 
ture or not, considered as an indication of their dynamical 
state, and we devised a substructure characterization of the 
haloes which is free of halo resolution effects. We then di- 
vide the haloes to those with and without a high level of 
substructure. 

Our results can be summarized in the following: 

(i) The well known relation between halo shape and halo 
mass has also an environmental dependence, albeit weak. 
Haloes found at small distances from their massive cluster 
show a systematic shift toward larger axial ratios (ie., they 
are more spherical) with respect to equal mass haloes at 
larger distances, a fact which is true for the whole range of 



halo masses studied. This result appear to be in disagree- 
ment with Hahn et al. (2006). 

(ii) The velocity dispersion of equal mass haloes shows a 
dependence on the environment. Haloes with substructure, 
near massive clusters, have a larger velocity dispersion with 
respect to equal mass haloes with no or low substructure 
index. This is probably due to the higher halo merging rate 
in high-density environments. On the contrary, haloes found 
further away do not exhibit this same behaviour: the veloc- 
ity dispersion- mass trend is the same independent of the 
presence or not of substructure. The velocity dispersion of 
high mass haloes does not seem to be affected by the en- 
vironment, nor that of any halo at a distance larger than 
~ 10ft~ 1 Mpc from its nearest massive cluster. 

(iii) The influence of environment is also reflected in the 
internal alignment of the velocity and density ellipsoid prin- 
cipal axes. Such an alignment is stronger for higher mass 
groups, probably due to the better definition of their shape 
given that these groups are more elongated than lower mass 
ones, while it is weaker near massive clusters, where the 
influence of the cluster and the high-density halo neighbor- 
hood is stronger. It is even weaker for haloes with a high 
level of substructure, which reflects the fact that during a 
merger the halo phase-space is significantly perturbed. 

(iv) Angular momentum and minor axes of haloes are 
roughly aligned, even more so for haloes with substructure. 
This relation does not seem to depend strongly on envi- 
ronment. However, one should keep in mind that the un- 
certainty of this measure is quite large, due to resolution 
effects. 

(v) On larger scales we detect alignments between the 
orientation of a halo and the direction to its nearest mas- 
sive cluster, which probably reflects the orientation of the 
filament within which they are embedded. The halo mi- 
nor/major axes appear perpendicular/parallel to the fila- 
ment, while the signal for both alignments is stronger for 
haloes near massive clusters and for haloes with no sub- 
structure. 

(vi) Overall we have found that the halo properties stud- 
ied in this work as a function of the distance to their nearest 
cluster, show a strong dependence on the amount of halo 
substructure. Since significant halo substructure is related 
to on-going or a recent merger, we could infer that the in- 
fluence of the close neighborhood of a halo, in the vicinity 
of massive clusters, is not less important than the influence 
of the cluster itself. 

There are at least two mechanisms involved in the evo- 
lution of the shape, alignment and velocity dispersion of 
haloes, namely the formation time (eg. Gottlober, Klypin & 
Kratsov 2001; Sheth & Tormen 2004) and the influence of 
the immediate environment. In this paper we were concerned 
with the latter aspect of the problem. We can summarize the 
interpretation of our results as follows. 

On the one hand, haloes forming in high density regions 
collapse earlier and they would on average have had more 
time to evolve and thus sphericalize, more so with respect to 
equal mass haloes forming in low-density regions. However, 
in the high-density environments an opposing factor is the 
overabundance of haloes which induce mergers and intra- 
halo interactions, which then disturb the virialized nature of 
these older haloes. Higher mass haloes evolve hierarchically 
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by the accretion of lower-mass haloes and thus are more 
elongated with respect to lower-mass haloes, which collapse 
and form earlier according to CDM models. 

The rising of the velocity dispersion of haloes (having 
significant substructure) near massive clusters, with respect 
to equal mass haloes with insignificant substructure or with 
those found in lower density regions, could be attributed to 
the higher halo merging rate present in the high-density en- 
vironment of massive clusters. The fact that the fraction of 
haloes with significant substructure is higher in high-density 
regions (see end of section 3), indeed reflects the more fre- 
quent halo mergers and interactions, which introduces also 
a bulk-flow (infall) velocity component in the halo velocity 
dispersion measure. 

The merging processes in high-density environments oc- 
cur along the anisotropic distribution of matter, which de- 
fines the large-scale filaments orientation. This is reflected 
in the alignments of the angular momenta of haloes, which 
are strongly influenced by the merging process, with the mi- 
nor axis of the halo and the alignment of the halo major axis 
with the orientation of the filament, defined by the direction 
between the halo and its nearest massive cluster. 



ACKNOWLEDGMENTS 

C.J. Ragone-Figueroa is supported by a CONICET fellow- 
ship, Argentina. This work has been partially supported by 
the European Commission's ALFA-II programme through 
its funding of the Latin-american European Network for 
Astrophysics and Cosmology (LENAC), the Consejo de In- 
vestigaciones Cientificas y Tecnicas de la Republica Ar- 
gentina (CONICET), the Secretaria de Ciencia y Tecnica 
de la Universidad Nacional de Cordoba (SeCyT), Agencia 
Nacional de Promotion Cientffica de la Republica Argentina 
and Agencia Cordoba Ciencia. 

REFERENCES 

Allgood B.F., Flores A., Primack J.R., Kravtsov A.V., Wechsler 

R.H, Faltenbachcr A., Bullock J.S., 2006, MNRAS, 367, 1781. 
Altay, G., Colberg, J.M., Croft, R.A.C., 2006, MNRAS, 370, 1422. 
Aragon-Calvo, M.A., van de Weygaert, R., Jones, B.J.T., van dcr 

Hulst, J.M., 2006, ApJ submitted, astro-ph/0610249. 
Avila-Reese V., Coh'n P., Gottlober S., Firmani C, Maulbetsch 

C, 2005, ApJ, 634, 51. 
Bailin J., Steinmetz M., 2005, ApJ, 627, 647B. 
Basilakos, S., Plionis, M., Yepes, G., Gottlober, S., Turchaninov, 

V., 2006, MNRAS, 365, 539. 
Bertschinger, E., 2001, ApJS, 137, 1. 

Bctt, P., Eke, V., Frcnk, C.S., Jenkins, A., Helly, J., Navarro, J., 
2006, MNRAS submitted, astro-ph/0608607. 

Boselli, A., Gavazzi, G., 2006, PASP in press, astro-ph/0601108. 

Bullock, J.S., 2002, in The Shapes of Galaxies and their Dark 
Matter Haloes, ed. Natarajan, P., World Scientific, p. 109. 

Colberg, J.M., Krughoff, S.K., Connolly, A. J., 2005, MNRAS, 
359, 272. 

Dressier A., 1980, ApJ, 236, 351. 

Dressier A. & Shcctman S.A., 1988, AJ, 95, 985. 

Dubinski J., 1992, ApJ, 401, 441. 

Einasto M., Einasto J., Miiller V., Hcinamaki P., Tucker D.L., 

2003, A&A, 401, 851. 
Einasto M., Suhhonenko I., Heinamaki P., Einasto J., Saar E., 

2005, A&A, 436, 17E. 



Evrard, A.E. et al., 2007, astro-ph/0702241. 

Faltenbacher, A., Gottlober, S., Kerscher, M., Miiller, V., 2002, 
A&A, 395, 1. 

Faltenbacher, A., Gottlober, S., Mathews, G., 2006, 
astro-ph/0609615. 

Goto, T., 2003, PhD Thesis, Univ. of Tokyo, astro-ph/0310196. 

Gottlober S., Turchaninov V., 2006, EAS Publications Series, Vol- 
ume 20, pp. 25-28, astro-ph/0511675. 

Gottlober S., Klypin, A., Kratsov, A.V., 2001, ApJ, 546, 223. 

Hahn, O., Porciani, C, Carollo, CM., Dekel, A., 2006, MNRAS 
submitted, astro-ph/0610280. 

Hopkins, P.F., Bahcall, N.A. & Bode, P., 2005, ApJ, 618, 1. 

Jing, Y.P & Suto, Y., 2002, ApJ, 574, 538. 

Kasun S.F. & Evrard A.E., 2005, ApJ, 629, 781K. 

Kazantzidis, S., Kratsov, A.V., Zcntner, A.R., Allgood, B., Nagai, 
D., Moore, B., 2004, ApJ, 611, L73. 

Knebe A. & Miiller V., 2000, A&A, 354J61K. 

Lee, J., 2006, ApJ, 644, L5. 

Lemson G. & Kauffmann G., 1999, MNRAS, 302, 111L. 

Ma, CP., & Bertschinger, E., 1995, ApJ, 455, 7. 

Maccio, A. V., Dutton, A. A., van den Bosch, F.C, Moore, 
B., Potter, D., Stadel, J., 2006, MNRAS submitted, 
astro-ph/0608157. 

Martinez H.J., Muriel H., 2006, MNRAS, 370, 1003M. 

Maulbetsch, C, Avila-Reese, V., Colin, P., Gottlober, S., 
Khalatyan, A., Steinmetz, M., 2006, ApJ, in press, 
astro-ph/0606360. 

Paz, D.J., Lambas, D.C, Padilla, N., Merchan, M., 2006, MN- 
RAS, 366, 1503. 

Pimbblet, K.A., 2005, MNRAS, 358, 256. 

Plionis, M., Ragone-Figueroa, C.J., Basilakos, S., 2006, in Groups 
of Galaxies in the Nearby Universe, eds. I.Saviane, V.Ivanov 
& J.Borissova. Springer- Verlag series "ESO Astrophysics 
Symposia", astro-ph/0604155. 

Plionis, M., 2004, in Outscirts of Galaxy Clusters: intense life in 
the suburbs, IAU Colloquium 195, cd. A. Diaferio. 

Plionis, M. & Basilakos, S., 2002, MNRAS, 329, L47. 

Ragone C.J., Merchan M., Muriel H., Zandivarez A. ,2004, MN- 
RAS, 350, 983. 

Schuecker, P., Boehringer, H., Reiprich, T. H., Feretti, L., 2001, 

A&A, 378, 408. 
Sheth R.K., Tormen C, 2004, MNRAS, 350 ,1385S. 
Springel V., 2005, MNRAS, 364, 1105S. 

Stoica, R., Martinez, V., Mateu, J., Saar, E., 2005, A&A, 434, 
423. 

Tormen, C, 1997, MNRAS, 290, 411. 

Vitvitska M., Klypin A., Kravtsov A., Wechsler R., Primack J., 

Bullock J., 2002, ApJ, 581, 799. 
van Haarlem M., & van de Weygaert, R., 1993, ApJ, 418, 544. 
West, M.J., 1994, MNRAS, 268, 79. 



© 0000 RAS, MNRAS 000, 000-000 



